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^ | Abstract. We provide a brief overview of some current experimental and theoretical 

. issues of heavy quark production. 
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The production of heavy quarks in high energy processes has become an increas- 
es) ■ ingly important subject of study both theoretically and experimentally. The theory 
of heavy quark production in perturbative Quantum Chromodynamics (pQCD) is 



INTRODUCTION 



more challenging than that of light parton (jet) production because of the new 
physics issues brought about by the additional heavy quark mass scale. The cor- 
rect theory must properly take into account the changing role of the heavy quark 
over the full kinematic range of the relevant process from the threshold region 
Q_i! (where the quark behaves like a typical "heavy particle" ) to the asymptotic region 

(where the same quark behaves effectively like a massless parton). With steadily 
improving experimental data on a variety of processes sensitive to the contribution 
of heavy quarks (including the direct measurement of heavy flavor production), this 
is a very rich field for studying different aspects of the QCD theory including the 
problems of multiple scales, summation of large logarithms, subtleties of renormal- 
ization, and higher order corrections. We shall briefly review a limited subset of 
these issues. 1 

THE FACTORIZATION THEOREM 

Perturbative calculations for heavy quark production are performed in the con- 
text of the factorization theorem expressed below in the commonly used form: 



f a ^ b (x, fi 2 ) ® <7 6 ^ C (Q 2 /// 2 , il4//i 2 , a s (fx)) + 0(A 2 QCD /Q 2 



*' Presented at 4th Workshop on Heavy Quarks at Fixed Target (HQ 98), Fermilab, Batavia, IL, 
10-12 Oct 1998. 

1 ) For a recent comprehensive review, see: Frixione, Mangano, Nason, and Ridolfi, Ref. [1] 
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FIGURE 1. Basic processes for DIS heavy quark production, a) 0(pP s ) flavor excitation: 
7 + c — > c; b) 0(a\) flavor creation: 7 + g — > c + c; c) 0(a|.) flavor excitation: 7 + c — ► c+ d) 
C(q!s) light-quark (q) fragmentation: (7 + g — »• q» + g) <g> (<? — ► c) . 



While the factorization was originally proven for massless quarks, [2] the theorem 
has recently been extended by Collins [3] to incorporate quarks of any mass, in- 
cluding "heavy quarks." (Note, we have explicitly retained the Mfj dependence in 
a.) It is important to note that the corrections to the factorization are only of 
order Aq CD /Q 2 , and not Mjj/Q 2 , even for the case of general quark masses. 

The factorization theorem can also be expressed as a composition of t-channel 
two particle irreducible (2PI) amplitudes: 2 
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Here, C represents the graph for a hard scattering, K represents the graph for 
a rung, T represents the graph that couples to the target, and Z represents a 
collinear projection operator. The first term in brackets roughly corresponds to 
the hard scattering coefficient function a, and the second term to the parton distri- 
bution function (PDF), /. Note that these two terms only communicate through a 
collinear projection operator, Z. Part of the effort in generalizing the factorization 
theorem for the case of massive quarks involves constructing the proper Z, and 
demonstrating that terms containing (1-Z) are power suppressed. However, once Z 
is determined, Eq.(2) yields an all-orders prescription for computing for both the 
hard scattering coefficient (a) and the parton distribution function (/). A calcu- 
lation using this formalism was first performed by ACOT [4] for the case of heavy 
quark production in deeply inelastic scattering, and we now examine this process 
in detail. 



HEAVY QUARK PRODUCTION IN DIS 

Several experimental groups [5] have studied the semi-inclusive deeply inelastic 
scattering (DIS) process for heavy-quark production £1 + N — > £2 + Q + X. New data 
from HERA investigates the DIS process in a very different kinematic range from 



2 ) I must necessarily leave out many details here; for a precise treatment, see Collins [3 
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FIGURE 2. The scaled cross section for DIS heavy quark production as a function of the quark 
mass rrifj- 

that available at fixed-target experiments. This perception has changed the way 
that we compute the semi-inclusive DIS heavy quark production. Traditionally, the 
heavy quark mass was treated as a large scale, and the number of active parton 
flavors was fixed to be the number of quarks lighter than the heavy quark. In this 
scheme, the perturbation expansion begins with the 0{ot\) heavy quark creation 
fusion process jg — » cc, (c/., Fig. lb). We refer to this approach as the Fixed 
Flavor Number (FFN) scheme since the number of flavors coming from parton 
distributions is fixed at three for charm production. 3 

More recently, a Variable Flavor Number (VFN) scheme (ACOT [4]) has been 
proposed which includes the heavy quark as an active parton flavor with non-zero 
heavy quark mass. In this case, the perturbation expansion begins with the O{a s ) 
heavy quark excitation process 7c — > c, (c/., Fig. la). The key advantages of this 
scheme are: [7] 

1. By incorporating the heavy quark into the parton framework, the composite 
scheme yields a result which is valid from threshold to asymptotic energies; in 
contrast, the FFN scheme contains unsubtracted mass singularities which will 
vitiate the perturbation expansion in the m c — ► or E — > 00 limit. 

2. Because the composite scheme resums the large logarithms appearing in the 
FFN scheme into the parton distribution functions, it includes the numerically 
dominant terms of the 0{a 2 s ) FFN scheme calculation in a 0(a l s ) calculation. 

In effect, the VFN scheme subsumes the FFN scheme. To illustrate this fact with a 
concrete calculation, in Fig. 2, we plot the cross section for "heavy" quark produc- 
tion as a function of the quark mass. 4 This figure clearly shows the three important 



3 ' The necessary diagrams have been computed to 0{aV) by Smith, van Neerven, and collabo- 
rators, c/., Ref. [6]. 

4 ) To be specific, we have computed single quark production for a photon exchange with x = 0.1, 
/i = Q = 10 GeV, and the cross section is in arbitrary units. 
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FIGURE 3. The contributions to DIS charged current inclusive F, 



charm 



vs. \i. For each separate 



contribution, the thick lines are the MS result (m s 
with m s = 0.5 GcV. From Kretzer, Ref. [8]. 



0), and the thin lines are the ACOT result 



kinematic regions. 1) In the massless region, where <C Q, the ACOT VFN re- 
sult reduces precisely to the massless MS result. 2) In the decoupling region, where 
n^H ^ Q, this "heavy quark" decouples and its contribution vanishes. 3) In the 
transition region, where m# ~ Q, this (not-so) "heavy quark" plays an important 
dynamic role. While the FFN scheme is appropriate only when m H ~ Q, we see 
that the VFN scheme is valid throughout the full kinematic range. 5 

This point is also illustrated in a calculation by Kretzer [8] (c/., Fig. 3) which 
shows the partial contributions to the charged current _P 2 c/iarm . 6 In this figure, each 
line is actually a pair of lines: the thin lines represents the result for _P 2 c/iarm using 
the ACOT scheme with m s = 0.5 GeV, and the thick lines regularize the strange 
quark with the massless MS prescription. (The charm mass is, of course, retained.) 
The fact that these two calculations match so closely (particularly in comparison 
to the pi- variation) indicates: 1) the ACOT scheme smoothly reduces to the desired 
massless MS limit as m H — * 0, and 2) for m H ^Aq C d we see that the quark mass 
no longer plays a dynamic role in the process and becomes purely a regulator. 



HEAVY QUARKS AND THE GLOBAL PDF ANALYSIS 

Recent precision data on F 2 and on F^ harm from HERA indicate that the charm 
contribution can rise to 25% of the total F 2 at small-x. These results clearly imply 
the need to perform new global analyses to account for the correct physics behind 

5 ) Buza et al, have determined the asymptotic form of the heavy quark coefficient functions 
which are then used to determine the threshold matching conditions between the three- and four- 
flavor shemes, Ref. [6]. Thorne and Roberts have a similar scheme with slightly different matching 
conditions, Ref. [9]. 

6 ) Kretzer and Schienbein have performed the first calculation of the 0(as) quark initiated 
process for general masses and general couplings, Ref. [8] . 
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FIGURE 4. Comparison of HI data in the small- a; region. From Lai and Tung, Ref. [10]. 



these measurements. Tung and Lai [10] have repeated the CTEQ4M global anal- 
ysis, [11] but this time implementing the heavy quark leptoproduction within the 
ACOT formalism to obtain a CTEQ4HQ set of PDF's. The deviation of CTEQ4HQ 
distributions from CTEQ4M are minimal, and are most noticeable at small-x; in- 
terestingly, the differences are larger for the light quarks than for the gluon and 
charm. 

The effect of these new PDF's and the comparison with data are shown in Fig. 4. 
The solid curves show the CTEQ4M distributions convoluted with massless ma- 
trix elements. The dashed curves show the CTEQ4M distributions convoluted 
with massive matrix element; while technically this is a mismatch of schemes, this 
comparison is useful to gauge the magnitude of the heavy quark effects, (which 
we observe are comparable to the experimental uncertainties). Finally, the dotted 
curves show the CTEQ4HQ distributions convoluted with massive matrix element. 
When a consistent scheme is used for both the matrix elements and the PDF's, the 
agreement with data is excellent. (This is as expected since this data was included 
in the fit.) It is interesting to note that overall x 2 fo r CTEQ4HQ (x 2 =1293) was 
slightly improved compared to the previous best fit CTEQ4M (x 2 =1320) for 1297 
data points. While this difference is small, we find it reassuring that the proper 
treatment of the heavy quark mass resulted in an improved fit; particularly when 
compared with a 4-flavor FFN fit (x 2 =1349) or a 3-flavor FFN fit (x 2 =1380). 

A recent re-analysis of the EMC data [12] concluded that there could be an 
intrinsic charm component in the proton of 0.86 ± 0.60%. It would be interesting 
to repeat this calculation in the context of a global analysis using the VFN ACOT 
formalism to see if more recent data favor an intrinsic charm component. 
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FIGURE 5. Comparison of F 2 from charged and neutral current DIS. From Seligman, et al, 
Rcf. [13]. 



HEAVY QUARKS AND EXTRACTION OF S(X) 

A topic closely related to DIS charm production is the extraction of the strange 
quark distribution. 7 In principle, we can extract s(x) by comparing DIS neutral 
and charged current data. To leading order, we have: 



F 2° 5 J\ 3(s + s)-(c + c) + ... 



F 2 CC 18 1 5 q + q 



(3) 



While the individual F 2 structure functions are measured precisely (c/., Fig. 5), [13] 
this approach is indirect in the sense that small uncertainties in the larger valence 
distributions will magnify the uncertainty on the extracted s(x). 

A direct method of obtaining s(x) is to use the neutrino induced di-muon process: 
u^N — > fi~cX with the subsequent decay c —> sfi + u^. Here, the di-muon signal is 
directly related to the charm production rate, which goes via the process W + s — >■ c 
at leading order. The method has the advantage that the signal from the s-quark 
is not a small effect beneath the valence process. 

A complete NLO experimental analysis was performed using the CCFR data 
set. [15] The recently collected data from the NuTeV experiment will provide an 
opportunity to extend the precision of these investigations still further. [16] Their 
high intensity sign-selected neutrino beam and the new calibration beam allows 
for large improvement in the systematic uncertainty while minimizing statistical 
errors. (See the paper by T. Adams, this meeting. [17]) 



7 ) For a comprehensive review, see Conrad, Shacvitz, and Bolton, Rcf. [14]. 
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FIGURE 6. The scaled cross section (nb-GeY 3 ) vs. p T for the a) LO-FFN and LO-VFN con- 
tributions, and b) TOT-FFN and TOT-VFN contributions. For each contribution, we choose 
fj, = [Mt/2,2Mt], with Mt = \J m 2 H + pj,, to gauge the /it-variation. From Ref. [18]. 

HADROPRODUTION OF HEAVY QUARKS 

We now turn to the hadropro duct ion of heavy quarks, and discuss how the 
method of ACOT [4,18] is used to provide a dynamic role for the heavy quark 
parton. We concentrate mostly on b-production at the Tevatron for defmiteness, 
and present typical results for b quark production. [1,19,20] (See the paper by 
A. Zieminski, this meeting. [21]) Fig. 6a shows the scaled differential cross sec- 
tion vs. pt for b production at 1800 GeV for the leading order (LO) calculations. 
The heavy creation (HC) process 8 (gg — > bb) represents the LO contribution to 
the fixed-flavor- number (FFN) scheme result. The heavy excitation (HE) process 
(gb —> gb) plus the HC term represents the LO contribution to the variable-flavor- 
number (VFN) scheme result. The pair of lines for each result shows the effect of 
varying \i. In a similar manner, Fig. 6b shows the total FFN and VFN results. 9 

Two interesting features are worth noting. 1) Examining Fig. 6a we observe 
the HE contribution is comparable to the HC one, in spite of the smaller 6-quark 
PDF compared to the gluon distribution. Closer examination reveals that two 
effects contribute to this unexpected result: a larger color factor and the presence 
of t-channel gluon exchange diagrams for the HE process, as compared to the HC 
process. 2) The LO-VFN (=HC+HE) contributions (Fig. 6a) (tree processes) give 
a reasonable approximation to the full cross section TOT-VFN (Fig. 6b); thus, the 
NLO-VFN correction is relatively small. This is in sharp contrast to the familiar 
FFN scheme where the TOT-FFN term is more than twice as large as the LO-FFN 
(=HC). This is, of course, an encouraging result, suggesting that the VFN scheme 

8 ) In this section we let g represent both gluons and light quarks, where applicable. Therefore, 
the HC process described as gg — > bb also includes qq — ► bb. 

9 > The formidable calculations of the NLO gg — > bb process were computed by Nason, Dawson, 
and Ellis (Ref. [22]), and also by Beenakker et at, (Ref. [23]). These calculations were implemented 
in a Monte Carlo framework (including correlations) by Mangano, Nason, and Ridolfi , (Ref. [24]). 
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FIGURE 7. F 2 vs. fi for DIS c-production. a) F 2 HE , F£ UB and the difference {F 2 HE - Fi UB ). 
The solid curves are for the mass-independent evolution scheme, and the dashed curves are for 
the mass-dependent evolution scheme, b) F™ T and F 2 HE — F 2 UB . The difference between 
the mass-independent evolution and mass-dependent evolution for F 2 is higher order and 
comparable or less than the [i- variation. From Ref. [27]. 

heavy quark parton picture represents an efficient way to organize the perturbative 
QCD series. 

In Fig. 6a, we also observe that while the TOT-VFN result provides minimal 
/x-variation for low pt, the improvement is decreased for large pt- This may be, in 
part, due to that fact that the TOT-VFN result shown here is missing the NLO-HE 
process gb — > ggb since this calculation, with masses retained, does not exist. In a 
separate effort, Cacciari and Greco [25] have used a NLO fragmentation formalism 
to resum the heavy quark contributions in the limit of large pt- This calculation 
effectively includes the massless limit of the gb — > ggb contribution (omitted above); 
the result is a decreased /i- variation in the large pt region. Recently, this calculation 
has been merged with the massive FFN calculation by Cacciari, Greco, and Nason, 
(Ref. [26]); the result is a calculation which matches the FFN calculation at lowpT, 
and takes advantage of the NLO fragmentation formalism in the high pt region. 

MASSIVE VS. MASSLESS EVOLUTION 

In a consistently formulated pQCD framework incorporating non-zero mass heavy 
quark partons, there is still the freedom to define parton distributions obeying 
either mass-independent or mass-dependent evolution equations. With properly 
matched hard cross-sections, different choices merely correspond to different fac- 
torization schemes, and they yield the same physical cross-sections. We demon- 
strate this principle in a concrete order a s calculation of the DIS charm structure 
function. [27] In Fig. 7 we display the separate contributions to F^ harm for both 
mass-independent and mass-dependent evolution. The matching properties are 
best examined by comparing the (scheme- dependent) heavy excitation F^ E and 
the subtraction F^ 113 contributions of Fig. 7a. 

We observe the following. 1) Within each scheme, F 2 HE and F 2 SC/S are well 
matched near threshold, cf., Fig. 7a. Above threshold, they begin to diverge, but 
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the difference (F^ E — F^ UB ), which contributes to -F 2 TOT , is insensitive to the 
different schemes. 2) It is precisely this matching of F BE and F 2 S ' C/B which ensures 
the scheme dependence of F^ 07 is properly of higher-order in a s , (cf., Fig. 7b). 

This matching is not accidental, but simply a result of using a consistent renor- 
malization scheme for both F% E and F2 UB '. To understand this we expand these 
terms near threshold (// ~ m#) where the mn/Q terms are relevant: 

VSUB = R f 9 ,P ® M^ce- = R f 9 /P ® ^ K R P { X ® + 

Here, the prescript i? specifies the renormalization scheme. From these relations, it 
is evident that F BE and i 7 ^ 5 * 75 will match to 0{a 2 s ) so long as a consistent choice or 
renormalization scheme R is made for the splitting kernels, R Pg^l c - This is the key 
mechanism that compensates the different effects of the mass-independent vs. mass- 
dependent evolution, and yields a a TOT which is identical up to higher-order terms. 
The lesson is clear: the choice of a mass-independent MS or a mass- dependent 
(non-MS) evolution is purely a choice of scheme, and becomes simply a matter of 
convenience- there is no physically new information gained from the mass- dependent 
evolution. 

CONCLUSIONS 

We have provided a brief overview of some current experimental and theoretical 
issues of heavy quark production. The wealth of recent heavy quark production 
data from both fixed-target and collider experiments will allow us to to extract 
a precision measurement of structure functions which can provide important con- 
straints for searches of new physics at the highest energy scales. As an important 
physical process involving the interplay of several large scales, heavy quark produc- 
tion poses a significant challenge for further development of QCD theory. 

We thank J.C. Collins, R.J. Scalise, and W.-K. Tung for valuable discussions, 
and the Fermilab Theory Group for their kind hospitality during the period in 
which part of this research was carried out. This work is supported by the U.S. 
Department of Energy, and the Lightner-Sams Foundation. 
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